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Abstract 



Results on two-particle angular correlations for charged particles emitted in pPb col- 
lisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV are presented. The 
analysis uses two million collisions collected with the CMS detector at the LHC. The 
correlations are studied over a broad range of pseudorapidity, tj, and full azimuth, 
(p, as a function of charged particle multiplicity and particle transverse momentum, 
px- In high-multiplicity events, a long-range (2 < |A?/| < 4), near-side (A^ 0) 
structure emerges in the two-particle Arj-Acp correlation functions. This is the first ob- 
servation of such correlations in proton-nucleus collisions, resembling the ridge-like 
correlations seen in high-multiplicity pp collisions at ^/s = 7 TeV and in A A collisions 
over a broad range of center-of-mass energies. The correlation strength exhibits a 
pronounced maximum in the range of px = 1-1.5 GeV/c and an approximately linear 
increase with charged particle multiplicity for high-multiplicity events. These ob- 
servations are qualitatively similar to those in pp collisions when selecting the same 
observed particle multiplicity, while the overall strength of the correlations is signifi- 
cantly larger in pPb collisions. 
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1 Introduction 

This Letter presents measurements of two-particle angular correlations in proton-lead (pPb) 
collisions at a nucleon-nucleon center-of-mass energy ^/s~^ = 5.02 TeV, performed with the 
Compact Muon Solenoid (CMS) detector at the Large Hadron Collider (LHC). Two-particle 
correlations in high-energy collisions provide valuable information for characterizing Quan- 
tum Chromodynamics and have been studied previously for a broad range of collision ener- 
gies in proton-proton (pp), proton-nucleus (pA), and nucleus-nucleus (AA) collisions. Such 
measurements can elucidate the underlying mechanism of particle production and possible 
collective effects resulting from the high particle densities accessible in these collisions. 

Studies of two-particle angular correlations are typically performed using two-dimensional 
Ari-A(p correlation functions, where A(p is the difference in azimuthal angle (p between the two 
particles and Arj is the difference in pseudorapidity rj =— ln(tan(0/2)). The polar angle 6 is 
defined relative to the counterclockwise beam. 

Of particular interest in studies of collective effects is the long-range (large \Ari\) structure of 
two-particle correlation functions, which is less susceptible to known sources of correlations 
such as resonance decays and fragmentation of energetic jets. Measurements in high-energy 
AA collisions have shown significant modifications of the long-range structure compared with 
minimum bias pp collisions IQI. Novel correlation structures extending over large Arj at |A^| ?a 
and I A^l ^ 2n/3 were observed in azimuthal correlations for intermediate particle transverse 
momenta, pj ~ 1-5 GeV/c pWTOj. In AA collisions, long-range correlations are interpreted as a 
consequence of the hydrodynamic flow of the produced strongly interacting medium lUTII and 
are usually characterized by the Fourier components of the azimuthal particle distributions 
[I2l. Of particular importance are the second and third Fourier components, called elliptic 
and triangular flow, as they most directly reflect the medium response to the initial collision 
geometry and its fluctuations ||T3| , and allow the study of fundamental transport properties of 
the medium using hydrodynamic models |(T4l - [T6l . 

In current pp and pA Monte Carlo (MC) event generators, the dominant sources of such long- 
range correlations are momentum conservation and away-side (A^ ^ n) jet correlations. Mea- 
surements in pp collisions at 7 TeV have revealed the emergence of long-range, near-side (A^ ^ 
0) correlations in a selection of collisions with very high final-state particle multiplicity lllTl . A 
large variety of theoretical models have been proposed to explain the origin of these so-called 
ridge-like correlations (see Ref. [18J for a recent review). The proposed mechanisms range from 
color connections in hard scattering processes and collective effects in the initial interaction of 
the protons to hydrodynamic effects in the high-density system possibly formed in these col- 
lisions. It is natural to search for the possible emergence of related features in pPb collisions, 
where a similar range of final-state multiplicities can be explored. The first comparison of pp 
and pPb measurements as a function of charged particle multiplicity and particle transverse 
momentum, presented in this Letter, should provide valuable information for understanding 
the origin of the long-range, near-side correlation signal seen in high-multiplicity pp collisions. 

2 Experimental Setup 

This analysis uses a pPb data set collected during a short run in September 2012. The beam 
energies were 4 TeV for protons and 1.58 TeV per nucleon for lead nuclei, resulting in a center- 
of-mass energy per nucleon pair of ^/s^ = 5.02 TeV. Due to the energy difference, the nucleon- 
nucleon center-of-mass in the pPb collisions is not at rest with respect to the laboratory frame. 
Since the higher energy proton beam traveled in the clockwise direction, i.e. atO = n, massless 
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3 Event and Track Selection 



particles emitted at rjcm = in the nucleon-nucleon center-of-mass frame will be detected at 
rj = —0.465 in the laboratory frame. 

A detailed description of the CMS experiment can be found in Ref. ||T9| . The main detector 
subsystem used for this analysis is the tracker, located in the 3.8 T field of the superconducting 
solenoid. It measures charged particles within the pseudorapidity range \rj\ < 2.5 and consists 
of 1440 silicon pixel and 15 148 silicon strip detector modules. It provides an impact parameter 
resolution of ~15 }im and a transverse momentum (pj) resolution of about 1.5% for lOOGeV/c 
particles. Also located inside the solenoid are the electromagnetic calorimeter (ECAL) and 
hadron calorimeter (HCAL). The ECAL consists of more than 75000 lead-tungstate crystals, 
arranged in a quasi-projective geometry and distributed in a barrel region (|?/| < 1.48) and two 
endcaps that extend up to \ri\ = 3.0. The HCAL barrel and endcaps are sampling calorimeters 
composed of brass and scintillator plates, covering < 3.0. Iron forward calorimeters (HE) 
with quartz fibers, read out by photomultipliers, extend the calorimeter coverage up to \ri\ = 
5.0. The detailed MC simulation of the CMS detector response is based on Geant4 pO|. 

3 Event and Track Selection 

The relatively low pPb collision frequency (about 200 Hz) provided by the LHC in this pilot 
run allowed the use of a track-based minimum bias trigger. Eor every pPb bunch crossing, the 
detector was read out and events were accepted if at least one track with pj > 0.4GeV/c was 
found in the pixel tracker. In the offline analysis, a coincidence of at least one HE calorimeter 
tower with more than 3 GeV of total energy on both the positive and negative sides of HE 
was required to select hadronic collisions. Events were also required to contain at least one 
reconstructed primary vertex within 15 cm of the nominal interaction point along the beam axis 
(Zvtx) and within 0.15 cm transverse distance to the beam trajectory. At least two reconstructed 
tracks were required to be associated with the primary vertex. Beam related background was 
suppressed by rejecting events with a high fraction of pixel clusters inconsistent with a single 
collision vertex [21]. A total of 2 million events passed all selection criteria, corresponding to 
an integrated luminosity of about 1 }ih^^, assuming a pPb interaction cross section of 2.1 barns. 

The angular correlation functions were obtained using the CMS highPurity f22] track selection. 
Additionally, a reconstructed track was only considered as a primary-track candidate if the sig- 
nificance of the separation along the beam axis, z, between the track and the primary vertex, 
dz/o'{dz), and the significance of the impact parameter relative to the primary vertex transverse 
to the beam, di/ cr{dj), were each less than 3. The relative uncertainty of the momentum mea- 
surement, a{pj) / Pj, was required to be less than 10%. To ensure high tracking efficiency and 
low fake rate, only tracks within \r]\ < 2.4 and with pj > 0.1 GeV/c were used. 

To match the analysis used for high-multiplicity pp collisions IITtII , the events were divided into 
classes of reconstructed track multiplicity, N°J^™, where primary tracks with \r]\ < 2.4 and 
Pt > 0.4GeV/c were counted. The fraction of events falling into each of the four multiplicity 
classes is listed in Table [ij The table also lists the average values of N°f^^^ and N^^^'^^'^, the 
event multiplicity of charged particles with | ?/ 1 < 2.4 and pj > OA GeV/ c corrected for detector 
acceptance and efficiency of the track reconstruction algorithm, as discussed in the following 
section. 
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Table 1: Fraction of the full event sample for each multiplicity class. The last two columns show 
the observed and corrected multiplicities, respectively, of charged particles with < 2.4 and 
Pj > 0.4GeV/c. Systematic uncertainties are given for the corrected multiplicities. 



Multiplicity class Fraction (w^f^'^J) ^JVj™'^'^ 

Minimum Bias 100.0 40.6 53.4±2.9 

j^offline < 35 50 4 17 1 23.5±1.3 

35 < N°f*™ < 90 41.9 56.3 75.6±4.1 

90 < N°f*"'' < 110 4.6 98.2 114.3±6.2 

jy^offline > 31 128.2 149.1±8.1 

4 Calculation of the Two-Particle Correlation Function 

The analysis of two-particle correlations was performed in classes of track multiplicity, 
following the procedure established in [7, 8|. For each track multiplicity class, "trigger" par- 
ticles are defined as charged particles originating from the primary vertex within a given pj 
range. The number of trigger particles in the event is denoted by Ntrig. In this analysis, particle 
pairs are formed by associating every trigger particle with the remaining charged primary par- 
ticles from the same pj interval as the trigger particle (a minimum of two particles is required 
in each pj bin from each event). The per-trigger-particle associated yield is defined as 

' """"" -B(0.0)x|<^. (1) 



Ntrig dA//dA^ ' B{Ar],A(p)' 

where Arj and are the differences in t] and (p of the pair. The signal distribution, S{Arj, A(p), 
is the per-trigger-particle yield of particle pairs from the same event, 

1 d^N^''""^ 

The mixed-event background distribution, used to account for random combinatorial back- 
ground and pair-acceptance effects, 

B(A?7,A<*) = " ' (3) 

Ntrig dA?/dA^ 

is constructed by pairing the trigger particles in each event with the associated particles from 
10 different random events in the same 2 cm wide Zytx range. The symbol N'"" denotes the 
number of pairs taken from the mixed event, while B(0, 0) represents the mixed-event asso- 
ciated yield for both particles of the pair going in approximately the same direction and thus 
having full pair acceptance (with a bin width of 0.3 in Arj and n/\6 in A(p). Therefore, the ra- 
tio B(0,0)/B(A?/, A^) is the pair-acceptance correction factor used to derive the corrected per- 
trigger-particle associated yield distribution. The signal and background distributions are first 
calculated for each event, and then averaged over all the events within the track multiplicity 
class. 

Each reconstructed track is weighted by the inverse of an efficiency factor, which accounts 
for the detector acceptance, the reconstruction efficiency, and the fraction of misreconstructed 
tracks. Detailed studies of tracking efficiencies using MC simulations and data-based methods 
can be found in |23|. The combined geometrical acceptance and efficiency for track reconstruc- 
tion exceeds 50% for pj ^ O.lGeV/c and |//| < 2.4. The efficiency is greater than 90% in the 
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5 Results 



|?/| < 1 region for pj > 0.6GeV/c. For the multiplicity range studied here, little or no depen- 
dence of the tracking efficiency on multiplicity is found and the rate of misreconstructed tracks 
remains at the 1-2% level. 

Simulations of pp, pPb and peripheral PbPb collisions using the PYTHIA, HIJING and HYDJET 
event generators, respectively, yield efficiency correction factors that vary due to the different 
kinematic and mass distributions for the particles produced in these generators. Applying 
the resulting correction factors from one of the generators to simulated data from one of the 
others gives associated yield distributions that agree within 5%. Systematic uncertainties due 
to track quality cuts are examined by loosening or tightening the track selections on dz/o'{dz) 
and dxy/o-{dxy) from 2 to 5. The associated yields are found to be insensitive to these track 
selections within 2%. 




Figure 1: 2-D two-particle correlation functions for 5.02 TeV pPb collisions for pairs of charged 
particles with 1 < pj < 3 GeV/ c. Results are shown (a) for low-multiplicity events (N°f ^'^'^ < 
35) and (b) for a high-multiplicity selection (N°^"^® > 110). The sharp near-side peaks from jet 
correlations have been truncated to better illustrate the structure outside that region. 



5 Results 

Figure [l] compares 2-D two-particle correlation functions for events with low (a) and high (b) 
multiplicity, for pairs of charged particles with 1 < pj < 3GeV/c. For the low-multiplicity 
selection (N"^™"^ < 35), the dominant features are the correlation peak near (Af/, A^) = (0,0) 
for pairs of particles originating from the same jet and the elongated structure at A^ ^ n for 
pairs of particles from back-to-back jets. To better illustrate the full correlation structure, the jet 
peak has been truncated. High-multiplicity events (N°^™'^ > 110) also show the same-side jet 
peak and back-to-back correlation structures. However, in addition, a pronounced "ridge"-like 
structure emerges at A^ extending to |Aj/| of at least 4 units. This observed structure is 
similar to that seen in high-multiplicity pp collision data at \/s = 7 TeV | fT7[ and in AA collisions 
over a wide range of energies ||3 HT0| . 

As a cross-check, correlation functions were also generated for tracks paired with ECAL pho- 
tons, which originate primarily from decays of tt'^s, and for pairs of ECAL photons. These 
distributions showed similar features as those seen in Fig. [l| in particular the ridge-like corre- 
lation for high multiplicity events. 

To investigate the long-range, near-side correlations in finer detail, and to provide a quanti- 
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Figure 2: Correlated yield obtained from the ZYAM procedure as a function of |A(^| averaged 
over 2 < | A?^ | < 4 in different px and multiplicity bins for 5.02 TeV pPb data (solid circles) and 
7 TeV pp data (open circles). The pi selection applies to both particles in each pair. Statistical 
uncertainties are smaller than the marker size. The subtracted ZYAM constant is listed in each 
panel. Also shown are pPb predictions for HIJING [24] (dashed curves) and a hydrodynamic 
model ||25j (solid curves shown for 1 < pj < 2GeV/c). 



6 



5 Results 



tative comparison to pp results, one-dimensional (1-D) distributions in are found by av- 
eraging the signal and background two-dimensional (2-D) distributions over 2 < \l^t]\ < 
4 13 muz]- In the presence of multiple sources of correlations, the yield for the correlation of in- 
terest is commonly estimated using an implementation of the zero-yield-at-minimum (ZYAM) 
method ||26| . A second-order polynomial is first fitted to the 1-D A^ correlation function in the 
region 0.1 < \A(p\ < 2. The minimum value of the polynomial, Czyam/ is then subtracted from 
the 1-D A^ correlation function as a constant background (containing no information about cor- 
relations) to shift its minimum to be at zero associated yield. The statistical uncertainty on the 
minimum level of obtained by the ZYAM procedure as well as the deviations found 

by varying the fit range in A^ give an absolute uncertainty of ±0.0015 on the associated yield, 
independent of multiplicity and pj. 

Figure |2] shows the results for pPb data (solid circles) for various selections in px and multiplic- 
ity Nj~, with px increasing from left to right and multiplicity increasing from top to bottom. 
The results for pp data at = 7TeV, obtained using the same procedure [IZJ, are also plotted 
(open circles). 

A clear evolution of the Acp correlation function as a function of both pj and is observed. 

For the lowest multiplicity selection in pp and pPb the correlation functions have a minimum at 
A(p = and a maximum at A^ = tt, reflecting the correlations from momentum conservation 
and the increasing contribution from back-to-back jet-like correlations at higher pj. Results 
from the HIJING [24J model (version 1.383), shown as dashed lines, qualitatively reproduce the 
shape of the correlation function for low N°^^'^. 

For multiplicities N"^"^*^ > 35, a second local maximum near |A^| emerges in the pPb 
data, corresponding to the near-side, long-range ridge-like structure. In pp data, this second 
maximum is clearly visible only for N°^™ > 90. For both pp and pPb collisions, this near- 
side correlated yield is largest in the 1 < px < 2GeV/c range and increases with increasing 
multiplicity. While the evolution of the correlation function is qualitatively similar in pp and 
pPb data, the absolute near-side correlated yield is significantly larger in the pPb case. 

In contrast to the data, the HijiNG calculations show a correlated yield of zero at A^ = for 
all multiplicity and px selections. The long-range, near-side enhancement is also absent in 
simulated pp collision events with the PYTHIA 1,27, 28 | event generator (version 6.4.24) and in 
simulated pPb collisions with the AMPT Il29ll model (version 1.25/2.25). 

Long-range correlations in pPb collisions have been quantitatively predicted in models assum- 
ing a collective hydrodynamic expansion Il25l . The correlation resulting from the predicted 
elliptic and triangular flow components for pPb collisions at y^s^ = 4.4 TeV are compared to 
the observed correlation in Fig. |2] for the 1 < px < 2GeV/c selection (solid line, second col- 
umn). The magnitudes for elliptic and triangular flow of V2 = 0.066 and Vj, = 0.037 correspond 
to those given in Ref. Il25ll for the highest multiplicity selection and the average value of px 
~ 1.4GeV/c found in the data. The same V2 and C3 coefficients were used for all multiplicity 
classes, showing the multiplicity dependence of the correlated yield assuming a constant flow 
effect. While this provides an indicative and useful illustration of the magnitude of the ob- 
served near-side enhancement, a detailed quantitative comparison of the model and data will 
need to include the additional non-hydrodynamical correlations from back-to-back jets, as well 
as the effects of momentum conservation, which suppress the correlation near A(p ^ relative 

to A(p ^ 71. 

The ridge-like structure in pPb collisions was also predicted to arise from initial state gluon 
correlations in the color-glass condensate framework ||30|. This model qualitatively predicts 
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the increase in the correlation strength for higher multiplicity pPb collisions, although it re- 
mains to be seen if the large associated yield seen in the highest multiplicity selection can be 
quantitatively reproduced in the calculation. 
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Figure 3: Associated yield for the near-side of the correlation function averaged over 2 < 
|A;y| < 4 and integrated over the region \A(p\ < 1.2 in 7 TeV pp collisions (open circles) and 
5.02 TeV pPb collisions (solid circles). Panel (a) shows the associated yield as a function of px 
for events with N°^™ > 110. In panel (b) the associated yield for 1 < pj < 2GeV/c is shown 
as a function of multiplicity The pj selection applies to both particles in each pair. The 

error bars correspond to statistical uncertainties, while the shaded areas denote the systematic 
uncertainties. 

The strength of the long-range, near-side correlations can be further quantified by integrating 
the correlated yield from Fig. |2] over |A^| < 1.2 using 12 classes of multiplicity. The resulting 
integrated "ridge yield", normalized by the width of the pj interval, is plotted as a function of 
particle pj and event multiplicity in Fig.|3]for pp (open circles) and pPb (solid circles) data. The 
error bars correspond to statistical uncertainties, while the shaded boxes indicate the systematic 
uncertainties. 

Figure |3|a) shows that the ridge yield for events with N"^'^*^ > 110 peaks in the region 1 < 
Pj < 2 GeV/ c for both collision systems. However, while the yield in pp collisions is consistent 
with zero for the 0.1 < pj < 1 GeV/c selection, it remains greater than zero in pPb data even 
for the 0.1 < px < 0-5 GeV/c range. For higher pj, the ridge yield in pp collisions is consistent 
with zero for px > 3 GeV/c, while the pPb results only approach zero at pj ~ 4-7GeV/c. 

The multiplicity dependence of the ridge yield for 1 < px < 2GeV/c particle pairs is shown in 
Fig. ^h). For low-multiplicity collisions, the ridge yield determined by the ZYAM procedure 
is consistent with zero, indicating that ridge-like correlations are absent or smaller than the 
negative correlations expected due to, e.g. momentum conservation. At higher multiplicity the 
ridge-like correlations emerge, with an approximately linear rise of the ridge yield observed 
for N°f ™ > 40, which corresponds to NX"''''''^ ^ 53. While the multiplicity dependence is 
qualitatively similar for pp and pPb collisions, a significantly larger yield per trigger particle is 
seen in pPb than in pp at a given multiplicity. 
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When interpreting the differences in the correlation structure between the two collision sys- 
tems, it is important to consider the relative contributions of different particle production 
mechanisms to the observed particle yields. While very high multiplicity pp collisions should 
mainly arise from rare multiple hard-scattering processes, the high-multiplicity pPb events 
should mostly result from particle production in multiple soft proton-nucleon scatterings. This 
will in particular affect the correlations due to back-to-back jet fragmentation in the A(p ^ n 
region. A simultaneous description of the measurements in pp and pPb should provide signif- 
icant constraints on models of the underlying physics processes. 

6 Conclusion 

The CMS detector at the LHC has been used to measure angular correlations between two 
charged particles with \r]\ < 2.4 in pPb collisions at ^^/s^ = 5.02 TeV. Azimufhal correlations 
for 2 < \Ar]\ < 4 in high-multiplicity pPb collisions exhibit a long-range structure at the near 
side (A^ 0). This ridge-like structure is qualitatively similar to that observed in pp collisions 
at ^/s = 7 TeV and in AA collisions over a broad range of center-of-mass energies. The effect is 
most evident in the intermediate transverse momentum range, 1 < pj < 1.5GeV/c. The near- 
side ridge yield obtained by the ZYAM procedure is foimd to be consistent with zero in the low- 
multiplicity region, with an approximately linear increase with multiplicity for N°f > 40 
(corresponding to js/^o^cted > 53^ While the multiplicity and pj dependences of the observed 
effect are similar to Ihose seen in pp data at = 7 TeV, the absolute ridge yield in pPb is 
significantly larger than in pp collisions of the same particle multiplicity. 
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